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We report on a new electromagnetic phenomenon that emerges in Mott insulators, i.e., materials that do not 
conduct electricity because of strong electronic Coulomb repulsion. The phenomenon manifests as antiferro- 
magnetic ordering due to orbital electric currents which are spontaneously generated from the coupling between 
spin currents and an external homogenous magnetic field. This novel spin-charge current effect provides the 
mechanism to detect the so far elusive spin currents by means of unpolarized neutron scattering, nuclear mag- 
netic resonance or muon spectroscopy. We illustrate this mechanism by solving a half-filled Hubbard model on 
a frustrated ladder, a simple but nontrivial case of strongly interacting electrons. 
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Insulators are generically characterized by the presence of 
a gap to charge-carrying excitations, and classified according 
to the origin of this gap. The charge gap of Mott insulators 
is driven by intra-atomic electron-electron Coulomb interac- 
tions, and the half-filled Hubbard Hamiltonian is the minimal 
model that describes its properties. Mott insulators exhibit a 
very broad spectrum of physical properties due to the corre- 
lated nature of their electronic localization. While electrons 
are completely localized in the Wannier orbitals of band insu- 
lators, Mott insulators always exhibit a partial electronic de- 
localization due to the finitude of the Coulomb repulsion, U, 
relative to the kinetic energy or hopping amplitude t. 

The combination of partial derealization with the 
fermionic nature of electrons leads to the well known anti- 
ferromagnetic exchange, J, between localized spins Sj. It 
was shown recently that this partial derealization also leads 
to electronic charge redistribution for a subclass of bond- 
ordered spin states (TJ. Even more surprisingly Mott insu- 
lators can also contain a distribution of orbital electric cur- 
rents that emerge for chiral spin orderings |Q] . These orbital 
currents are directly associated with the notion of scalar spin 
chirality HO 

Xjkl = SjXS k -Si. (1) 

The main purpose of this Letter is to demonstrate that a uni- 
form magnetic field can induce orbital antiferromagnetism, 
i.e., antiferromagnetic ordering due to orbital electric cur- 
rents, by the stabilization of staggered spin currents and the 
subsequent interplay between these currents and the applied 
field. For this purpose we will consider a Hubbard Hamil- 
tonian on a zig-zag chain (see Fig[T| with nearest-neighbor 
(next-nearest-neighbor) hopping amplitudes t\ (£2). The ap- 
plication of a magnetic field B increases the uniform magne- 
tization along the field direction z. A spontaneous breaking of 
the remaining U(l) symmetry of global spin rotation along the 
z-axis is prevented by quantum fluctuations (Mermin- Wagner 
theorem [4Q. This implies the absence of usual magnetic or- 
dering in the x-y plane: (Sj) = with 77 = x, y. We will 
show however that the so called "chiral spin orderings", both 
vector chirality ((Sj x S J+1 ) ^ 0) and scalar spin chirality 



((Xijk) 7^ 0), can be induced by the field. The field-induced 
vector chiral ordering is a remnant of the classical helical or- 
der (Sj) ^ that is obtained in the classical (large spin, 
S — > 00) limit |5|. On the other hand, the scalar spin chi- 
rality results from the spontaneous vector spin chirality and 
the field induced magnetization, (xjkl) — (Sj x Sj.) ■ (Sj). 
As we will show below, this connection between scalar and 
vector spin chiralities implies a relation between electric and 
spin currents, i.e., a spin-charge current effect. 

Chiral phases (nonzero scalar or vector spin chirality) in 
quantum spin chains were predicted a long time ago and have 
been studied for many years. As noted by Villain [6|, the chi- 
ral ordering must survive at finite temperature (T) in weakly 
coupled chains, without the usual helical long-range ordering, 
since the chiral correlation length is much longer that the spin 
correlation length. In other words, there must be a window 
of temperatures where chiral order takes place in absence of 
helical order. There are experimental indications of the ex- 
istence of such phase in the quasi-one dimensional organic 
magnet Gd(hfac)3NITiPr [5 7|. However, the absence of ex- 
ternal physical fields that couple to non-uniform chiral order- 
ings poses a challenge for measuring these phases when there 
is no helical ordering [8]. By showing that the field-induced 
vector chirality is accompanied by a staggered ordering of 
orbital magnetic moments (orbital electric currents), we pro- 
vide a direct way of detecting this exotic phase by means of 
unpolarized neutron scattering, nuclear magnetic resonance 
(NMR) or muon spectroscopy ( /i-SR). 

We start by considering a half-filled Hubbard model on the 
zigzag ladder with L sites (0 < j < L — 1) depicted in Fig{T] 
and in the presence of magnetic field B = Bz 

j,v,a j 

where v = 1, 2, er = ±1, ct(c JCT ) creates(annihilates) an 

electron of spin a at site j, rij a = c i j cr Cj a ., rij = J2 a an< ^ 

Sj = Ea^ c ja a !p c jp/2 are trie spin- 1/2 operators with a v 
the Pauli matrices. We assume periodic boundary conditions 
(PBC), i.e., L = 0. In the large U/t v limit, the low-energy 



spectrum of H is described by an effective Heisenberg spin- 
1/2 Hamiltonian, 

H = J2MSj-Sj+v- \)~ B Y, S h ( 3 ) 

because the electrons are localized near the lattice sites. H 
is obtained by projecting the original H into the low-energy 
subspace S. This is true in general for any physical quan- 
tity, A, whose effective low-energy operator, A, is a func- 
tion of the spin operators Sj. The expression for A is ob- 
tained by a canonical transformation that follows from stan- 
dard degenerate perturbation theory. The exchange constants 
are J v = 4±l/U > 0. 
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FIG. 1: Zigzag ladder. The arrows indicate the circulation of the spin 
and the electrical orbital currents that emerge in the ground state for 
J i ^> Ji and \B\ < |B sat |, with _B sat the saturation field. The small 
circles indicate the orientation of the (staggered) magnetic moments 
generated by the electric currents. 

We are interested in two physical quantities that will be 
particularly relevant for this work. These are the charge and 
the spin (z-component) current densities in a given bond (jl) 
(with I = j + v and tji = t„) 

cr 

r n = *X](4<^-4s'a)^- f ^ ( 4 ) 

cr 

where f 3 ; = (r; — r 3 )/|r; — rj\, and Yj is the vector position 
of site j. Both are Noether currents for the charge and spin 
conservation laws associated with the corresponding U(l) and 
SU(2) global symmetries of H. The charge current has a non- 
zero low-energy effective operator, whenever the site j 
belongs to a loop that is closed by an odd number of hopping 
terms (TJ. Since the shortest possible loop is a triangle, 1^ 
is 0(t n+1 /U n ) with n odd and n > 3. Therefore, the low- 
est order contribution to the effective charge current density 
operator is [ 1 1 

h ~ ~Y v o l 1^ tj2 Xna- ( 5 ) 

Equation |5]l establishes a direct correspondence between the 
scalar chiral spin ordering, (Xijk) ^ 0, and the presence of 
electric orbital currents. In other words, most of the scalar 
spin orderings will be accompanied by the emergence of or- 
bital currents. The effective spin current density operator is 

% = y4p with k ji = s j x s ' ' ( 6 ) 



This simple expression shows that the spin current is directly 
associated to the vector spin chirality. Since is even under 
a particle-hole transformation (tji — > —tji, Sj — > — Sj), the 
prefactor on the right hand side of Eq. (|6]l can only contain 
even powers of the hopping amplitudes. 

The effective current density operators introduced above 
are necessary to characterize the ground state correlations 
when the magnetic field approaches its saturation value B sa t, 
and t-z > t\. As we will show below, this ground state exhibits 
long range order of spin and charge currents that are roughly 
proportional to each other when J2 ^> J\, and \B\ < |-B S at|- 
To understand the origin of this instability, it is convenient to 
write H in terms of new fermionic degrees of freedom fj by 
means of a Jordan- Wigner transformation 

S t = f] K v Sj = K, J f v n^Sf + i (7) 

where fij = fjfj and Kj = Yik<j(^ ~ 2^fe) is the nonlocal 
operator that realizes the change in exchange statistics. 
The new expression of H can be written as 

H = Hx + Hz-BY^nu (8) 

1 

Hi = y J2(f}fl+1 + f} + ifl ) + Jl ]>>^+l " »0 

1 1 

where nf = n z ll+l = i(f}fj, +1 - fi+Ji)- Interestingly, we 
remark that the first term of H 2 is an explicit ferromagnetic in- 
teraction between spin currents on adjacent bonds. However, 
a state with net nearest-neighbors spin currents, (nf ) ^ 0, can 
only appear when Ji 7^ 0, i.e., for a finite coupling between 
upper and lower chains (see Fig{T|. 

To develop some intuition about the role played by J\, it is 
convenient to rewrite H in momentum space 

H = J2(tk - M)4«fe + v pi a l+P a l-p a k a q ^ ( 9 ) 

k kqp 

with a\ = Y^Zq e lfcj 7j, £fc = Ji cos k + J 2 cos 2k, fi = 
B + J1 + J2, and v pq = 2e p — 8J2 cos (p + 2q). The first term 
is the non-interacting part of H, while the second contains the 
density-density interactions of Hi and H2 that lead to the first 
contribution, 2e p , plus a correlated second-nearest-neighbor 
hopping that is contained in the first term of Hi, and leads to 
the second contribution in v pq . 

For J2 > Ji/4, the fermion dispersion tk has two de- 
generate minima at k = ±Q, with cosQ = —J1/AJ2, that 
correspond to opposite values of n\. The saturation field 
B sat = Ji + J 2 — €q (— i? S at) corresponds to the critical 
value of the chemical potential at which the fermion density 
p = (fij) becomes equal to one (zero). From now on we 
will assume that the spin system is close to full polarization. 
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Since the physics is independent of the sign of B (B — > — B 
under a time reversal transformation), we will choose the neg- 
ative sign, B > —B aat , for the rest of the analysis. For the 
equivalent fermionic problem, this condition is equivalent to 
the dilute density limit, p -C 1 or /i > eq. 

For Ji = 0, the symmetric and the anti-symmetric sec- 
tors generated by the fermionic operators a\ ± a\ , w are per- 
fectly decoupled. In this case the spin currents are quenched 
since nf can only connect states on different sectors (it is 
odd under the transformation k — > k + n). The situation 
changes dramatically for J1/J2 <C 1. At low energies, we 
can only have fermions at near the two degenerate minima: 
k ~ ±Q. Since Q is close to tt/2, the fermions <2j_q gen- 
erate spin currents of nearly maximal amplitude and opposite 
signs: 2~Zz K ? a l|0) = sm k a. J. 1 0) . While the non-interacting 
part of H favors a ground state in which both minima are 
equally populated (Pauli exclusion principle), the second term 
of v pq leads to a nearest-neighbor attraction between fermions 
with the same chirality (in the same minimum) and a nearest- 
neighbor repulsion between fermions with opposite vector 
chiralities (different minima). This is also expected from the 
first term of Hi and it implies the possibility of a chiral in- 
stability similar to the one originally proposed by Nersesyan 
and coauthors [9| for XY zig-zag spin chains (see also ifTUll ). 
In fact, the chiral phase was recently found for J 2 3> J\ and 
\B\ ~ _B sat [5] by using the mean-field decoupling of the 
bosonized version of H that was originally introduced in Ref. 
[9 1 for the anisotropic case. On the other hand, a straightfor- 
ward mean-field approximation to Eq. Q overestimates the 
stability of the chiral phase (it gives a wrong density depen- 
dence for the energy of the disordered state). Since mean-field 
approximations to interacting quasi-one-dimensional systems 
are always subjected to scrutiny, it is decisive to have a numer- 
ical confirmation of this phase. By using density-matrix renor- 
malization group (DMRG), McCulloch and coauthors ifTTIl re- 
cently found the chiral phase for J 2 = J\, while Okunishi 
|[T2l reported a phase diagram that confirms the existence of a 
chiral phase for J 2 ^> Ji and \B\ ~ J5 sat . 

We will now derive an important physical consequence of 
a chiral spin phase that was overlooked before. To simplify 
notation we introduce the definitions: xi = Xu+11+2 an d 
= 1/ with 7 = c, s. The combination of a field induced 
magnetization, m z — (S*), and a net vector spin chirality 
(nf ) 7^ leads to a non-zero mean value of the scalar spin 
chirality 
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Figure [T] displays the circulation of electrical orbital currents 
(arrows) that are expected for the chiral-ordered ground state 
of the ziz-zag chain (J 2 ^> J\ and \B\ < |B S at|) according to 
Eq. 
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Remarkably, the application of a uniform magnetic 
field induces orbital antiferromagnetism via the (dominant) 
Zeeman coupling to the spin moments. The small circles in 
Fig. [T] indicate the orientation of the staggered orbital mo- 
ments that these electric currents generate. 

To test the validity of Eq. ( [TT] ) we study the ground state 
properties of both H and H by means of the DMRG method 
lfT3l . We use PBCs 1141 to eliminate spurious oscillations in 
the correlation functions due to boundary effects. The incom- 
mensurate nature of the ordered ground state makes the nu- 
merical calculation quite challenging |[T5l . In contrast to Ref. 
lfT2ll we do not include any infinitesimal bias field and com- 
pute chiral-chiral correlators (instead of the order parameter 
(ki)) to establish chiral long-range order. In the case of the 
Hubbard Hamiltonian H , we solve a chain of L — 64 sites for 
t 2 /t 1 = 1.6, U/h = 20,24, and 30. B is chosen such that 
2m z — 0.875. The charge-current, C c (r) = (Iflf +r ), and 
the spin-current, C s (r) = (Iflj 1 ^), correlation functions are 
computed directly in the DMRG ground-state wave function. 
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FIG. 2: Size dependence of the two-point vector spin chirality corre- 
lator C s nm {r) for J2/J1 = 2.56 and 2m z = 2{Sf ) = 0.875. 



Xi ^ m z [«f + «f +1 - k, z u+2 ] . (10) 

This simple expression shows that, to a good a approximation, 
scalar and vector chiralities are proportional to each other, and 
the uniform magnetization m z is the proportionality constant. 
This relationship attains a clear physical meaning when we 
look at Eqs. <|5j and ([6j: the field induced vector chiral order 
contains orbital electric currents that are proportional to the 



For the Heisenberg model H, we study chains of lengths 
L — 48, 64, and 128 sites. The ratio between the exchange 
constants is determined by the ratio between the hopping am- 
plitudes J 2 /Jx = t\jt\ = 2.56 and 2m z = 0.875. In this 
case we compute the vector spin chirality two-point corre- 
lation functions: C^ m (r) = (Kii +n ni +r i +r+m ). Figure [2] 
shows Cfi(r) for different system sizes. The finite-size scal- 
ing clearly indicates the presence of long-range vector spin 
chiral order. 
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According to Eq. (|6j, the effective spin-current two-point 
correlation function is C s (r) — [8^-] 2 C'f 1 (r). At long dis- 
tances, the effective charge-current correlation function can 
be obtained approximately by using Eq. 



10) 



C c (r) ~ &(xiXi+r) * 4&m 2 z [C s n (r) - Cf 2 (r) + 



Cf a (r), 



with <5 = 2304tjtl/U 4 . Figures |Ja), (b) and (c) com- 
pares the spin-current two-point correlators computed with the 
Heisenberg (H) and the Hubbard (H) models for three differ- 
ent values of U/h = 20, 24 and 30. 

As expected the results for the two models show better 
agreement as U increases. According to Eq. ( [12} , the long- 
range ordered spin-currents (see Fig. [2} must lead to long- 
range ordered orbital electric currents. Figures |3ja), (b) and 
(c) show a comparison between C c (r) — (If If, r ) computed 
with the Hubbard Hamiltonian (H) for the same three differ- 
ent values of U, and the right-hand side of Eq. ( 12 1 computed 



with the Heisenberg model H. The good agreement between 
both curves confirms that the spin and the electric currents 
are approximately proportional to each other, with the propor- 
tionality constant linear in m z . In other words, the ordering 
of spin currents is accompanied by ordering of orbital electric 
currents for non-zero m z as follows from Eq. dTOb . 
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FIG. 3: (a), (b), (c) Charge-current correlation function C c (r) and 
(d), (e), (f) spin-current correlation function C s (r) in the case of the 
Hubbard model H (dashed lines) and the Heisenberg model H (solid 
lines). The results are shown for L = 64, tz/h = 1.6 (J2/J1 = 
2.56), (sf) = 0.4375, and U/h = 20, 24, and 30. 

The resulting staggered configuration of orbital magnetic 
moments implies that the ground state under consideration is 
a field-induced orbital antiferromagnet. We note that the or- 
bital moments are located at the centers of the triangles (see 
Fig. [TJ, while the spin moments are obviously located at the 
corners (or sites of the lattice). The values of the orbital cur- 
rents that we obtained for U/h = 20 (see Fig. ^d)) are of 
order 0.01 t\. The corresponding orbital magnetic moment is 
a few percent of a Bohr magneton for h ~ leV. In contrast to 
the spin currents (/ef), these orbital magnetic moments can be 



measured with different experimental techniques such as neu- 
tron scattering, NMR or /i-SR. This provides a simple way of 
detecting this exotic spin ordering in real materials. Moreover, 
notice that the magnitude of the magnetic moment would be 
much larger if the chiral phase remains stable in the interme- 
diate coupling regime U > t v . 

To induce an electric current out of an spontaneously gen- 
erated spin current (or vice-versa) one needs an externally ap- 
plied magnetic field (or any other field that is odd under a 
time reversal transformation) since both currents have oppo- 
site parity under time reversal. This symmetry consideration 
has to be complemented by a microscopic mechanism that de- 
termines the magnitude of the spin-charge current conversion. 
The possibility of having orbital electric currents in frustrated 
Mott insulators [ 1 1 is the key for finding and quantifying such 
mechanism. This spin-charge current effect cannot be associ- 
ated with a magnetoelectric response because the spin current 
does not couple directly to any electric or magnetic fields. On 
the other hand, the induced orbital electric currents can lead to 
small magnetoelectric effects. For instance, the application of 
an electric potential between the upper and lower chains (see 
Fig. [TJ will lead to a net orbital magnetization because the 
magnetic moments in the lower and upper triangles will not 
be compensated any longer. 

This work was carried out under the auspices of the Na- 
tional Nuclear Security Administration of the U.S. Depart- 
ment of Energy at Los Alamos National Laboratory under 
Contract No. DE-AC52-06NA25396. 
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